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Abstract 

Gravity takes care of both inflation and subsequent reheating in Starobinsky's R 2 - 
model. The latter is due to inflaton gravitation decays dominated by scalar particle 
production. It is tempting to suggest that dark matter particles are also produced in 
this process. Since free scalars being too hot cannot serve as viable dark matter [3], 
we further study the issue considering two options: scalars with non- minimal coupling 
to gravity and superheavy scalars generated at inflationary stage. We found that the 
first option allows for viable warm or cold dark matter if scalar mass exceeds 1.1 MeV. 
The second option implies supercold dark matter with particle mass 10 16 GeV, which 
production is saturated at the end of inflation when inflaton-dependent scalar mass 
rapidly changes and violates adiabaticity. Similar result holds for superheavy fermion 
dark matter. 

Cosmology of the homogeneous and isotropic Universe in Starobinsky's R 2 -mode\ [1] has 
an inflationary stage. This stage is realized as a large-field slow-roll inflation driven by 
scalaron field, that is a scalar degree of freedom coming from gravity sector of the theory. 
Subsequent decay of the scalarons into ordinary matter dominated by scalar particles pro- 
duction reheats the Universe. It is natural to speculate that dark matter particles could be 
produced in the same way. 

Indeed, it is gravity — the mostly universal force — that provides interaction for the scalaron 
with all (matter) fields. The strength of coupling to a given field is determined by the strength 
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of conformal symmetry violation. This leads to the fact that mass terms, being conformally 
noninvariant, become modulated by the scalaron field. Moreover, scalars have extra coupling 
due to their conformally-noninvariant kinetic term. Consequently, the scalaron decays are 
dominated by scalar particles production [1, 2], so that free scalars produced in an appro- 
priate amount would be too hot and could not serve as the dark matter [3]. At the same 
time, free heavy fermions can do the job. However, scalar dark matter implies only one new 
degree of freedom to be added to the SM, and therefore is the most economic solution. 

There are two obvious ways to get rid of the problem with scalar dark matter and avoid 
over-energetic scalar relics. First, one can partially restore the conformal symmetry by adding 
nonminimal coupling of scalar to gravity. By changing the value of corresponding coupling 
constant £ one can control the strength of conformal symmetry violation and govern the 
scalaron decay rate into dark matter particles. Then heavier and hence colder scalar relics 
become viable dark matter. The second option is to consider superheavy scalar particles 
which kinematically can not be produced from scalaron decays at postinflationary stage. 
However, some amount of them can be generated nonperturbatively in the time- dependent 
scalaron background. 

In this work we consider both options. We obtain the mass of the scalar dark matter 
candidate as a function of nonminimal coupling to gravity £. Depending on the value of £ 
dark matter produced at postinflationary epoch can be light or heavy, the lightest viable free 
scalars are of 1.1 MeV and form warm dark matter. Such particular candidate is interesting 
in the context of emerging problems of the standard CDM cosmology on small scales, such as 
missing satellites, galactic matter density profiles and angular momentum of spiral galaxies, 
for review see e.g. [4]. 

Considering the second option we find that the production of superheavy particles is 
saturated at the end of inflation, when the scalar mass modulated by the scalaron field 
changes rapidly. A nontrivial result is that the mass dependence of the number density of 
created particles is a power law. Thereby in this model it is possible to produce superheavy 
relics with mass greatly exceeding the scalaron mass, which form supercold dark matter. 

We start by recalling some relevant for the present study facts about i? 2 -inflationary 
model. The gravitational sector of this model is described by the following Lagrangian 

where R is scalar curvature, /x is dimensionful parameter, and we introduce the reduced 
Planck mass M P = M P i/\/&k = 2.4 x 10 18 GeV. Variation of the action (1) with respect to 
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9fiu yields the fourth order differential equation. One can argue (using general coordinate 
invariance) that the degrees of freedom of the field g^ u can be split into massless spin-2 field 
g^y and massive scalar field <fi which obey the second order equations of motion [5]. This 
splitting is achieved explicitly in the action written in the Einstein frame, which one can 
come to by conformal transformation 

&«/-»■ 9ixv = X9ixv , X = exp ^y^2/3(f)/Mp^ . (2) 
For the action (1) we have then 

M l 5 . ^ , 3 // 2 M- 



S EF = Lf=$d*x 



2 R+ 2^ d ^-^T JL \ 1 -m 



(3) 



everywhere g^ v stands for the metric tensor and R is scalar curvature for the metric g^ v . 
One observed that in the Einstein frame the action splits into the Einstein-Hilbert action 
f° r 9yLv and ordinary action for the scalar field <fi with a specific potential. The scalar mode 
in the action (1) was named scalaron [1] and we will use this name for the field in what 
follows. 

The scalaron potential in (3) is exponentially flat at super-Planckian field values and 
provides the slow-roll inflation. Normalization of the amplitude of scalar perturbations, 
generated during this stage, to the observed CMB anisotropy requires [6] 

\i = 1.3 x HT 5 M P . 

The spectral index and parameters of the generated tensor perturbations (for recent discus- 
sion see [7]) are consistent with observational constraints [8]. 

When the slow-roll conditions get violated, inflation terminates and the inflaton <fi starts 
to oscillate rapidly, with frequency equal to scalaron mass fi. This drives the Universe 
expansion like at matter-dominated stage. The intermediate stage naturally ends up with 
inflaton decays into ordinary particles. Here we briefly discuss this process, see [1, 2, 3] for 
details. 

The major role in the scalaron decay is played by scalars. Below we consider scalar (p 
described by the following Lagrangian in the Jordan frame 



Sf = j^d'x Q gTd^ip - \ my + ) , (4) 



where the last term represents nonminimal coupling to gravity with dimensionless parameter 
|£| < 1. Note that with |£| < 1 we do not introduce a new scale. Below we assume that £ is 
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a free small parameter, but let us single out two special points: minimal coupling with £ = 
and conformal coupling corresponding to £ = 1/6. The first value leads to a free scalar field 
theory while the second one ensures conformal invariance of the action (4) without mass 
term. The mass term explicitly breaks conformal invariance and one can expect (small) 
deviation of £ from the conformal value due to quantum corrections. Quantum corrections 
can also generate small but nonzero value of £ for a scalar field minimally coupled to gravity. 
Therefore one can expect that the value of coupling £ should be within some intervals covering 
£ = or £ = 1/6 rather than be exactly equal to £ = or £ = 1/6, respectively. Below 
we are interested in only positive values of nonminimal coupling. Negative £ may lead to a 
nonadiabatic evolution and thereby an amplification of scalar perturbations on superhorizon 
scales during inflation [9]. That in turn may change predictions for scalaron mass /i and 
requires detailed study beyond the scope of this paper. 
After conformal transformation (2) and field rescaling 



V? ->■ <P = X V 2 <p 



the action (4) takes the form 



~ ( 5 ) 

One can find a conformal transformation involving tp and leading to another Einstein frame, 
where the nonminimal coupling disappears (see, e.g. [9]). This frame is more convenient to 
study the case of two-field inflation, where the dynamics of scalar field ip is also relevant at 
the inflationary epoch. However, below we are interested in a more simple situation, where 
only scalaron is responsible for the inflation, and adopt the metric (2). 

Equation (5) describes the scalaron interaction with free scalars. It doesn't vanish at 
£ = and provides scalaron decay into a pair of scalars. The second term in Eq. (5) 
yields on nontrivial cosmological background a contribution to the scalar mass. However, at 
reheating one has \R\ <C fi 2 , and the scalaron decay into scalars of mass m v < fi/2 remains 
kinematically allowed 1 . Its width reads [2] 



r € =(l-« + ^) V 192,M2 • ^ 



1 At the inflationary stage the situation is different, as we discuss later on. 
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Massless scalar field with £ = 1/6 is conformally invariant and does not interact with the 
scalaron, which is illustrated by nullifying decay rate (6) at £ — > 1/6, m v — > 0. The mass 
term in Eq. (5) explicitly breaks conformal invariance and the scalaron decay rate in this 
case is proportional to m*, see Eq. (6). Hence, if the mass of the conformally coupled scalars 
is very small, m v <C /i, the scalaron decays very slowly. 

The value of the coupling constant, which is natural and guarantees the successful re- 
heating through scalaron decays, is £ = 0. In case of the Standard Model with four scalar 
degrees of freedom (in the Higgs sector) minimally coupled to gravity scalaron rapidly decays 
to scalars thereby reheating the Universe upto the temperature [1, 2, 3] 



- reh 



3.1xl0 9 GeV. (7) 



Now let us consider a new scalar field, which does not interact with SM particles or this 
interaction is very weak, so these particles never equilibrate in the primordial plasma. If 
stable at cosmological time scale, such scalar is a good candidate to be dark matter. Previ- 
ously [3] we have shown that minimally coupled scalars are too hot to be dark matter. Here 
we rehabilitate light scalars as a possible dark matter candidate by introducing nonminimal 
interaction with gravity. 

As we discussed above, light scalars are eventually produced from scalaron decays. In 
order to describe this process we will treat the scalaron field configuration as a condensate 
of nonrelativistic particles with number density n^t) and energy density p^{t) = fin^t). At 
some moment the number density of scalar particles produced from scalaron decays is 

To get the number density at reheating time t re h one should take into account the volume 
expansion, 

dn(t reh ) = dn(t m )-^K = 2 %*(Q4^\dt* ■ (8) 

a 6 {t r eh) P ^ {treh) 

Now let us look at the physical momentum of these particles. At production time t* it equals 



P* — yV 2 / 4 — m|. Then it gets redshifted and at reheating time one has 

Preh = P* °f *\ ■ (9) 

O-ytreh) 

Using this formula one can find the momentum difference at reheating time for the particles 
produced from t* to t* + dt*, 

dpreh = PrehH(t*)dt* . (10) 
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a. 




Figure 1: The distribution function for the scalars produced from scalaron decays. 



On the other hand one can introduce the distribution function f(p re h) of produced particles 
as dn(t re h) = fipreh) d 3 p re h/C27i) 3 . Substituting Eqs. (8), (10) into this formula, we get 



fip 



reh 



4?r 2 pS*) a3 (**) 



where production time i* should be expressed in terms of momentum by Eq. (9). At t re h 
one rewrites Eq. (11) as 

T 5 M 2 P H(Q 



fip 



reh ) 



P 



Pi 



(12) 



As one can see from Eq. (11), for calculation of the spectrum one needs to know the time 
dependencies of the scale factor, Hubble parameter and scalaron energy density. In order 
to obtain a more accurate expressions, we use numerical calculations. Our result for the 
spectrum of produced scalars is presented in Fig. 1. 

For the particles number density one obtains then 

n v ~ 2.5 ^M 2 P H reh , (13) 
P 

where numerical factor in front of r.h.s. of Eq. (13) is the result of numerical integration of 
the spectra in Fig. 1 over momenta. For the average momentum one gets 



p reh = 0.85 . 



(14) 
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Note that p re h ^> T re h, so the particles are very hot in comparison with the plasma. 

In order to find the present abundance of the dark matter one uses the ratio of the entropy 
to the number density. At the reheating time we have 

n(t reh ) ' r e M P ' v ; 

where g* is the effective number of degrees of freedom in the plasma, see e.g. [10]. The ratio 
(15) remains intact during the subsequent hot stages of the Universe expansion including 
the present epoch. From the requirement that the relative contribution of the dark matter 
to the present energy density p c is [8] Qdm = 0.223 one gets the following relation between 
the particle mass and nonminimal coupling £ which enters Eq. (6), 

m ^0.2 nDMpcn ^ reh . (16) 
s T(.M P 

Finally let us find the lower and the upper bounds on the scalar mass. The lower bound 
corresponds to the case of warm dark matter. These particles should have the average 
velocity v™ ax ~ 10~ 3 at the epoch of equilibrium between radiation and matter densities. 
On the other hand, the average velocity at that time can be found from Eq. (14). That leads 
to the particle mass 

0-42 I Q* : eq\ 3 T e q 



,max \ n IT 
u eq \ y* / ± i 



— ^-1-lMeV, (17) 



- reh 



where at the equilibrium T eq 0.76 eV, g* m ~ 3.9 and = 106.75 [10]. Then Eq. (16) 
implies £ = 1/6 ± 0.018 for the viable dark matter candidate. 

The upper bound can be imposed if we take into account the gravitational particle pro- 
duction from vacuum fluctuations in the expanding Universe. It has been shown [11, 12] that 
the particles conformally coupled to gravity and heavier than m^^^ ~ 10 9 GeV overdose 
the Universe. Possible dark matter masses (16) together with the lower and upper bounds 
are shown in Fig 3. 

The aforesaid consideration is valid while the occupation number for produced particles 
remains small, f{p re h) ^ 1- Otherwise, coherent effects become important and the corre- 
sponding occupation number undergoes exponential growth. Let us check the maximum 
occupation number, which in our case is reached by the particles produced at the end of 
inflation as one can see from Eq. (12). Using the relation (16) one finds that the maximal 
occupation number exceeds unity only for the very light dark matter particles with 

m v < 1 keV . (18) 
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These values of are smaller than (17), hence the coherent effects become significant for 
the irrelevant hot dark matter only, and our result (16) remains intact. 

We proceed to the case, when dark matter particle is heavier than the half-scalaron. 
Perturbative production of such particles is kinematically forbidden. However, some amount 
of them is eventually produced in a time-dependent scalaron background. In order to de- 
scribe this process one can use the method based on the Bogoliubov's transformation coeffi- 
cients [13]. Here we summarize the basic formalism which we employed in our analysis. For 
more details see, e.g Ref. [14]. 

We start by canonically quantizing the action (5) in curved cosmological background 
with the external classical scalaron field. For this purpose it's convenient to choose the 
FLRW metric in the form ds 2 = a 2 (i])(dri 2 — dx 2 ), where a(rj) is the scale factor and rj is the 
conformal time defined as drj = dt/a. After rescaling the field variable, 

<P = s / a (v) , 

the equation of motion following from the action (5) becomes 

id 2 d 2 1 2 2 (I \ ( a" <j)' 2 VQa 2 dVU)\} , ^ 

{ a? - W + x a < V<3 ) \~a + Jl} + ~M~p a<6 j j s ' r '' ' • < 19 > 

where prime means the derivative with respect to conformal time and V((j>) is the scalaron 
potential, see Eq. (3). Here we have exploited the classical equation of motion to express the 
second conformal time derivative of the scalaron field. Solution of Eq. (19) can be written 
in the following form 

*(V, 5) = — j d 3 p {a p s p ( V )e~^ + als;( V )e^) , (20) 

where d p and d p are annihilation and creation operators, and the mode function s p (r]) obeys 
the oscillatory equation 

S ; + u 2 (r ] )s p = (21) 

with time-dependent frequency 

2 , . 2 1 2 2 f l A ( a" (j)' 2 VQa 2 dVU)\ , , 

w (i?) =p + x v ~ \6 ~ / y~a + Ti}, + ~m~ f J ' < 22 » 

While to varies adiabatically, | u'\/u 2 1, solution (20) corresponds to a particle-like field 
configuration (note that adiabatic field evolution is a necessary condition for introducing 
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the very idea of particles). At this stage the particle creation does not occur. The number 
density of dark matter particles is given by 

^=(2^)3/ ^| 2 > ( 23 ) 

where (3 P is the Bogoliubov's transformation coefficient, which relates initial and final particle- 
like states of the quantum field. It can be expressed in terms of the mode functions as 



|2 



|s'| 2 + W 2 |c| 2 1 



(24) 

Assuming initially no dark matter in the Universe we impose the vacuum initial conditions 

s p — > 1 / V2u , s' p ^-—iujs p , at 77 — > — 00 . (25) 

As mentioned above, the particle number is conserved during adiabatic evolution. On 
the contrary, highly varying frequency gives rise to particle production and the higher the 
variation rate the larger the number of produced particles. In our case the maximal violation 
of the adiabaticity occurs right at the end of inflation due to the scalaron-dependent mass 
term. Particle creation in the case of large m v can be especially sensitive to scalaron dynamics 
at this stage. Therefore we resort to numerical methods in order to find exact evolutions of 
the scalaron field and the scale factor, and further solve numerically Eqs. (21), (22) for the 
mode function. Because of the limited capability of numerical calculations we impose the 
initial conditions (25) not at the past infinity but at some finite value of rj. We found that the 
results become practically independent of the choice of the initial time corresponding to more 
than 7 — 8 e-foldings till the end of inflation. It was also learned from numerical simulations 
that the number of particles in comoving volume produced at the end of inflation remains 
almost constant during subsequent evolution. Indeed, the number of produced particles is 
very small, hence no any Bose-enhancement is expected during the subsequent evolution. On 
the other hand, both inflaton velocity and the Hubble parameter slow down as the Universe 
evolves at the post inflationary matter-dominated stage. As one can see from (22), this leads 
to strengthening of adiabaticity and very low rate of particle production after the end of 
inflation. 

A convenient quantity for subsequent computation is the ratio of the energy density of 
produced scalar particles to the total energy density. This ratio remains almost constant 
while the Universe expands due to oscillating scalaron. Numerical results for this ratio in 
the case of minimally coupled scalars, £ = 0, are presented in Fig. 2 (small circles). It turns 
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Figure 2: The ratio of the energy density of produced scalar particles minimally coupled to 
gravity to the total energy density Small circles represent numerical results, while the line 
corresponds to numerical fit with the formula (26). 

out that numerical data are well described by the relation 



which is plotted in Fig. 2. Here the value of a for a given £ we obtained by fitting numerical 
data with the formula (26). We found that a varies from 10 -14 to 10~ 16 , when £ varies 
from to 1/4. Note, that the resulting number density of produced particles exhibits power 
law mass-dependence. This differs from the cases of superheavy particle production by 
pure gravity (see, e.g. [11, 12]) or by combined effects of gravity and inflaton field with 
polynomial coupling to the scalars (see, e.g. [15]), where the mass dependence of particles 
number density is exponential as mass exceeds expansion rate. In our case at inflationary 
stage the effective mass of scalars determined by the inflaton field is very small: i.e. they are 
not superheavy. Number density of produced at this stage scalars are exponentially decreases 
due to inflation. Scalars become superheavy at <p < Mp when inflation terminates. The 
scalar mass depends exponentially on the value of scalaron field (see Eq. (5)). Consequently 
the particle production is saturated at the end of inflation, when scalaron-dependent mass 
rapidly changes from <C H to > H violating adiabaticity, and exponential decrease 




(26) 



10 



in the number density due to the Universe expansion does not occur. From this chain of 
reasonings we conclude that it is exponentially-changing mass term which is responsible for 
the power law mass dependence (26). It would be worth, however, to find more rigorous 
arguments (analytic calculations) to support the numerical result (26). 

Making use of Eq. (26) and the requirement that produced scalar particles should explain 
dark matter phenomenon, one gets the relation between the dark matter mass and the 
nonminimal coupling constant, 



The lighter particles are overproduced in this model in accordance with Eqs. (26), (27). 
However, for very light scalars, < /i/2, Eq. (26) is not applicable. We find numerically 
in this case that the scalar production rate actually goes down with decreasing m v so that 
overproduction of particles with mass of two or three orders of magnitude lighter than the 
scalaron does not occur. However in this region of masses the dominant process is the 
perturbative production at reheating stage given by Eq. (6). 

Our final results for possible masses of scalar dark matter given by Eqs. (16), (27) are 
shown in Fig. 3 by solid thick lines. The top of the figure corresponding to superheavy dark 
matter particles with > /i/2 is plotted in linear scale. Here the ^-dependence of dark 
matter mass is almost linear as one can see from the figure. The filled area above the curve 
represents cosmologically allowed masses for any free scalars in this model, which do not 
overdose the Universe. 

The case of light dark matter which masses are given by Eq. (16) is presented in loga- 
rithmic scale. The zoom shows the fine structure of the curve near conformal point £ = 1/6. 
The filled region under the curve represents allowed masses of free scalar particles in the 
model. Allowed region is bounded from above by dashed line corresponding to scalar mass 
of 10 9 GeV. This bound appears when we consider all possible mechanisms of particle pro- 
duction which are relevant for this masses. Indeed, from Fig. 3 one concludes that the cosmo- 
logically allowed free scalars with m v < 10 12 GeV are only those which are almost conformal. 
On the other hand, the scalar particles conformally coupled to gravity are produced in ex- 
panding Universe from vacuum fluctuations and overdose the Universe if m v > 10 9 GeV [12]. 
The mass range 10 12 GeV < < /i/2 is excluded, if we take into account particle produc- 
tion during inflation, see the discussion below Eq. (27). Note, that limits obtained in [12] 
are not applicable to the model with nonminimal coupling outside the conformal window 
£ w 1/6 and scalar mass approaching /i/2 from below, see Fig. 3, which we do not study 




(27) 
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Figure 3: Mass of the scalar dark matter candidates as a function of nonminimal coupling 
constant £ (solid thick lines). The filled area represents cosmologically allowed masses of 
free scalars. At the bottom part of the plot the allowed region is located between horizontal 
lines corresponded to ~ 10 9 GeV and ~ 1.1 MeV. The zoom shows the fine structure 
of the curve near the conformal point £ = 1/6. 



in this work. Also, light free scalars of <C IkeV may be forbidden for £ in particular 
regions, where the exponential amplification of scalar production is expected due to coherent 
effects, see Eq. (18) and discussion therein. 

The special choices of dark matter parameters are marked by numbers in Fig. 3. The 
points with label "1" correspond to warm dark matter with mass given by Eq. (17). The 
point "2" represents conformally coupled scalars with ~ 2.8 x 10 15 GeV. The free dark 
matter scalars have ~ 1.3 x 10 16 GeV, which is marked by "3" in Fig. 3. 

Note in passing, that free superheavy fermions are also produced during inflation. We 
calculate their number density as we did for scalars. In order to be dark matter these 
fermions should have mass ~ 3.1 x 10 15 GeV, which is similar to the case of conformally 
coupled scalars presented above. It is worth to emphasize, that superheavy relics produced 



12 



at the end of inflation, if unstable, cannot help to reheat the Universe earlier. As we have 
seen from Fig. 2, these particles take away a very small fraction of the energy. One can also 
evaluate their number density to entropy ratio. For example, for the fermions of 10 12 GeV at 
reheating we have n^/s ~ a ~ 10~ 15 . This value is much smaller than t]b- Therefore, seesaw 
sterile neutrinos produced by this mechanism cannot help with leptogenesis and one must 
exploit the perturbative production of seesaw sterile neutrinos by the scalaron decay [3] , to 
explain the neutrino oscillations and the baryon asymmetry of the Universe. 

To conclude, we have shown that in i? 2 -inflation the scalar particles can form dark matter, 
which may be warm, cold or supercold depending on scalar mass and the value of nonminimal 
coupling constant to gravity (see Fig. 3). 

We end up with several comments. First of all, note that introduction of selfinteracting 
terms for the scalars can change our predictions for dark matter masses. Indeed, for the scat- 
tering processes which preserve the particles number, the average momentum is conserved 
and our predictions remain unchanged. However, if selfinteraction is strong enough so that 
scattering to multi-particle final states becomes rapid, the particle number density increases 
while average momentum goes down. In order to keep constant the relative contribution of 
the dark matter to the present energy density, one should decrease the mass of the candi- 
date for a given value of £. This decreasing factor for the mass is equal to the decreasing 
factor for the average momentum so that the average velocity at the epoch of equilibrium 
between radiation and matter densities would be the same. This in particular means that 
the boundary values of nonminimal coupling constant (marked by "1" in Fig. 3) will refer to 
the warm dark matter case as before. 

Second, superheavy dark matter, if slowly decaying, can show up in ultra-high energy 
cosmic rays. The cosmologically long but finite lifetime can be explained by a nonperturba- 
tive mechanisms of their decays, such as in the instanton scenarios [16], or involving quantum 
gravity (wormhole) effects [17]. 
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